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We discuss a fully relativistic Landau Fermi liquid theory based on the Quantum Hadro- 
Dynamics (QHD) effective field picture of Nuclear Matter (NM). From the linearized kinetic 
equations we get the dispersion relations of the propagating collective modes. We focus our 
attention on the dynamical effects of the interplay between scalar and vector channel con- 
tributions. A beautiful "mirror" structure in the form of the dynamical response in the 
isoscalar/isovector degree of freedom is revealed, with a complete parallelism in the role re- 
spectively played by the compressibility and the symmetry energy. All that strongly supports 
the introduction of an explicit coupling to the scalar-isovector channel of the nucleon-nucleon 
interaction. 

In particular we study the influence of this coupling (to a 5-meson-like effective field) 
on the collective response of asymmetric nuclear matter (ANM). Interesting contributions 
are found on the propagation of isovector-like modes at normal density and on an expected 
smooth transition to isoscalar-like oscillations at high baryon density. 

Important "chemical" effects on the neutron-proton structure of the mode are shown. 
For dilute ANM we have the isospin distillation mechanism of the unstable isoscalar-like 
oscillations, while at high baryon density we predict an almost pure neutron wave structure 
of the propagating sounds. 
PACS: 24.10.Cn, 24.10.Jv, 21.30.Fe, 24.30.Cz 



I. INTRODUCTION 



The QHD effective field model represents a very succesfull attempt to describe, in a fully consistent 
relativistic picture, equilibrium and dynamical properties of nuclear systems at the hadronic level |j]-0|. Very 
nice results have been obtained for the nuclear structure of finite nuclei @ |], for the NM Equation of State 
and liquid-gas phase transitions and for the dynamics of nuclear collisions . Relativistic Random- 
Phase- Approximation (RRPA) theories have been developed to study the nuclear collective response [fU3|-[f5[ . 

In this paper we present a relativistic linear response theory with the aim of a trasparent connection 
between the collective dynamics and the coupling to various channels of the nucleon-nucleon interaction. In 
particular we will focus our attention on the dynamical response of asymmetric nuclear matter since one of 
the main points of our discussion is the relevance of the coupling to a scalar isovector channel, the virtual 
<5[ao(980)] meson, not considered in the usual dynamical studies. Another point of interest is the dynamical 
treatment of the Fock terms, neglected in the usual Relativistic Mean Field (RMF) scheme of the papers 
cited before. 

A relativistic extension of the Landau linear response theory of Fermi Liquids has been considered before 
just starting from the relativistic form of the Landau parameters |lt|-|l8| . We will show that the full dispersion 
relations obtained from the relativistic kinetic equations present some interesting corrections that cannot be 
neglected. 

The main physics results are: 

• The important effect of a <5-meson coupling on the isovector collective mode at saturation baryon 
density. This is of interest for the relativistic study of the Giant Dipole Resonance in heavy finite 
nuclei. It is important to note that the inclusion of Fock terms is acting in the same direction. 

• The presence of noticeable "chemical effects" in the propagating collective oscillations, i.e. the charge 
symmetry of the "waves" is quite different from the asymmetry of the initial equilibrium matter. The 
effect is opposite for the unstable modes present at low densities, more proton rich and leading to the 
isospin distillation effect, and for the stable propagating sounds at high baryon density which appear 
mostly like "pure neutron waves" . 
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• A stimulating "mirror" structure of the isoscalar and isovector linear response, with the restoring forces 
given by the potential part respectively of the compressibility and the symmetry energy. The interplay 
between the scalar and vector meson effective fields in the dynamics is very similar for the two degrees 
of freedom, as already observed for static properties |l9| . The conclusion is that a consistent relativistic 
effective field model has to include on the same footing isoscalar and isovector meson fields, both scalar 
and vector. 

Our results around normal density can be used as general guidelines in predicting the behaviour of volume 
collective modes in finite /3-unstable nuclei. Similar study for asymmetric NM have been performed in 
Ref. [^0| using Skyrme-like interactions. Apart the difference in the used interactions, in particular for the 
symmetry terms, we will see similar results and interesting new relativistic effects. 

In Sect. II we derive the kinetic equations in the general case of non-linear self-interacting terms, including 
the Fock corrections. We discuss the inclusion of the <5-meson channel, also on the model parameters. In 
Sect. Ill we present the relativistic linear response equations. In Sect. IV we have a general discussion on the 
formal structure of the dispersion relations, the role played by the scalar/vector mesons and the comparison 
to non-relativistic cases. Results for isovector(-/ifce) collective modes are presented in Sect.V, in particular for 
the asymmetry and baryon density effects. The isoscalar(-/zfce) response is analysed in Sect. VI. Conclusions 
and outlooks can be found in Sect. VII. 



II. KINETIC EQUATIONS FROM A QHD EFFECTIVE THEORY 

We start from the QHD effective field picture of the hadronic phase of nuclear matter . In order to 
include the main dynamical degrees of freedom of the system we will consider the nucleons coupled to the 
isoscalar scalar a and vector u) mesons and to the isovector scalar 5 and vector p mesons. 

The Lagrangian density for this model, including non-linear isoscalar/scalar cr-terms plfl, is given by: 



C = ^(ifl" - g^ - g p B» ■ r) - (M - g a <f> - g s r ■ 6)ty 
~ m 2 s 4> 2 ) - ^4? - -^f - \w p „W^ + \mlV v V v 



l -(d,8 ■ d»8 m 2 6 2 ) -G pv ■ + \ m 2 p B v ■ B" (1) 

where W» v {x) = d*V v (x) - d^V^x) and G» v (x) = d^B v (x) - d v B»(x) . 

Here tp(x) is the nucleon fermionic field, (f>(x) and V v {x) represent neutral scalar and vector boson fields, 
respectively. S(x) and B v (x) are the charged scalar and vector fields and T denotes the isospin matrices . 

From the Lagrangian Eq. (Q) with the Euler procedure a set of coupled equations of motion for the meson 
and nucleon fields can be derived. The basic approximation in nuclear matter applications consists in 
neglecting all the terms containing derivatives of the meson fields, with respect to the mass contributions. 
Then the meson fields are simply connected to the operators of the nucleon scalar and current densities by 
the following equations: 

+ A$ 2 + 5$ 3 =$(x)il>(x) =pS (2) 

V"{x) = f^(xh^(x) = f u % , 
B"(z) = f p i;(x)rr^(x) , 
S(x) = M(x)tiP(x) (3) 

where $ = g^tj), f a = {g a /m a ) 2 , A = a/gl, B = b/g*, f w = (fc/m u ) 2 , f p = {g p /2m p ) 2 , f s = (g 5 /ms) 2 . 

For the nucleon fields we get a Dirac-like equation. Indeed after substituting Eqs.(H]|) for the meson 
field operators, we obtain an equation which contains only nucleon field operators. All the equations can be 
consistently solved in a Mean Field Approximation (RMF), where most applications have been performed, 
in particular in the Hartree scheme [||J^] . 

The inclusion of Fock terms is conceptually important [ p2|j23| since it automatically leads to contributions 
to various channels, also in absence of explicit coupling terms. We will discuss this point later. A thorough 
study of the Fock contributions in a QHD approach with non-linear self-interacting terms has been recently 
performed |2if| , in particular for asymmetric matter p5|| . 

The present approximation implies that retardation and finite range effects in the exchange of mesons 
between nucleons are neglected. Nevertheless, thanks to the small Compton wave-lengths of the mesons cr, 
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ui,p and J, the assumptions expressed by Eqs. (§,0) are quite reasonable. For light mesons such as pions this 
approximation is not justified. However in this case a perturbative expansion in the pion-nucleon coupling 
constant seems to be reasonable p2[ . Moreover it has been shown that the inclusion of pions does not change 
qualitatively the description of nuclear matter around normal conditions ]22] ] . 

We remark that the kinetic approach discussed here is fully consistent with the previous approximation. We 
are concerned with a semiclassical description of nuclear dynamics, so that the nuclear medium is supposed 
to be in states for which the nucleon scalar and current densities are smooth functions of the space-time 
coordinates. 

Within a mean field picture of the QHD model we focus our analysis on a description of the many-body 
nuclear system in terms of one-body dynamics. This is enough for the scope of the paper. Correlation effects 
can be effectively included at the level of coupling constants, as noted in the discussion of the results. 

We will perform the many-body calculations in the quantum phase-space introducing the Wigner transform 
of the one-body density matrix for the fermion field |?6|,^7j ■ 

The one-particle Wigner function is defined as: 

[F{x,p)U = 7^4 / d*Re-^ R (: ^{x + ^ a (x f ) :> , 

where a and (3 are double indices for spin and isospin. The brackets denote statistical averaging and the 
colons denote normal ordering. The Wigner function is a matrix in spin and isospin spaces; in the case of 
asymmetric NM it is useful to decompose it into neutron and proton components. Following the treatment 
of the Fock terms in non-linear QHD introduced in Ref. |^,^5|, we obtain for the Wigner function the 
following kinetic equation: 



fMxW ± f P 33»(x)'f - f*Ps(x)T fsPssix)} F (l) (x,p) = 0, i = n,p (4) 
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where A = d x ■ d p , with d x acting only on the first term of the products. Here ps3 = psp — psn and 
33[i{x) = jfi( x ) ~ iu(^) are the isovector scalar density and the isovector baryon current, respectively. We 
have defined the kinetic momentum and effective masses, as: 

M* = M — Up s (x) ± f sPS3 (x) , (5) 
with the effective coupling functions given by: 



/;4-S- ^rP^fm + k + b. - b 



p s 8dp s (x) 2p s y 'dp 2 s (x) 2 JW 2' 
_ 1 d*(x) 5 3 3 



J6 



- 1 d£(x) 1 _1 9 
U 8d Ps (x) + 2 IuJ 2 Jp+ ° JS 



~ _ 1 d*(x) 1 3 1 
Ip 8dp s (x) + 4 /w+ 4 /p 8 JS ' 

(6) 

where 8 TrF 2 (x) = pg + j^j 7 ' + p| 3 + J3^J 3m - We remind that we are dealing with a transport equation so 
the currents and densities, in general, are varying functions of the space-time, at variance with the case of 
nuclear matter at equilibrium. 

The expression of Eq.(||) for the effective mass, embodies an isospin contribution from Fock terms also 
without a direct inclusion of the 6 meson in the Lagrangian. The usual RMF approximation (Hartree level) 
is covered by the Hartree-Fock results, one has has only to change the coupling functions fi(i = a,u>,p,5), 
Eqs.(^|), with the coupling constants /j. 



Equilibrium properties: the nuclear Equation of State 

In the following we will study the collective modes. In order to analyze the results it is essential to relate 
them to the equation of state (EOS), that we will briefly discuss in the following. In particular for the 
collective modes in asymmetric nuclear matter it is important the behaviour of the symmetry energy E sym . 
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The energy density and pressure for symmetric and asymmetric nuclear matter and the n,p effective 

2 

masses can be self-consistently calculated just in terms of the four boson coupling constants, fi = (-%), 
i = a, to, p, S, and the two parameters of the a self-interacting terms, A = and B = see ref. f2-i| , |2~5|| . 

The isoscalar meson parameters are fixed from symmetric nuclear matter properties at T = 0: satu- 
ration density pa = 0.16/m -3 , binding energy E/A = — lQMeV, nucleon effective mass M* — 0.7Mn 
(Mjv = 939MeV) and incompressibility Ky = 240MeV at p . The fitted f a ,f u ,A, B parameters are re- 
ported in Table I. They have quite standard values for these minimal non-linear RMF models. Set I and 
Set II correspond to the best parameters within a non-linear Hartree calculation, respectively with the 
p — (Set I, NLH + p) and with the p + S (Set II, NLH + (p + 5)) couplings in the isovector channel (see 
the discussion in ref. jl9]]). NLHF stands for the non- linear Hartree- Fock scheme described before. 



Table I. Parameter sets. 



parameter 


Set I 


Set II 


NLHF 


NL3 


U (fm z ) 


11.27 


same 


9.15 


15.73 


U (/m a ) 


6.48 


same 


3.22 


10.53 


f P (/m a ) 


1.0 


2.8 


1.9 


1.34 


/* ifm") 


0.00 


2.0 


1.4 


0.00 


A (fm-') 


0.022 


same 


0.098 


-0.01 


B 


-0.0039 


same 


-0.021 


-0.003 



In the table we report also the NL3 parametrization, widely used in nuclear structure calculations pB[ . 
We remind that the iVL3-saturation properties for symmetric matter are chosen as po — 0.148/m -3 , M* 
ii.fi \ / v . K v = 271.8MeV. The symmetry parameter is a 4 = 37.4Mey. 

The symmetry energy in ANM is defined from the expansion of the energy per nucleon E(pB, a) in terms 
of the asymmetry parameter a defined as 

q, = _P^± _ PBn ~ Pb p - N ~ Z 
Pb Pb A 

We have 

E(p B ,a) = e SE^l = E(p B ) + E sym ( PB )a 2 + 0(a 4 ) + ... (7) 

PB 



and so in general 



l d 2 E(p B ,a) l _ 1 dh 
Esym = 2 drf _ 2 PB dpl~ 3 lpB3=0 (8) 



In the Hartree case an explicit expression for the symmetry energy can be easily derived p9| . [L9| 

E sym (PB) = - ^ + ^fpPB - 2 fs E* F 2 [l + f s A(k F: M*)} = E ^ m + EV °y m (9) 



where kp is the nucleon Fermi momentum corresponding to p B , E* F = y (kp + M* 2 ) and M* is the effective 
nucleon mass in symmetric NM, M* — Mm — g a 4>. 
The integral 

A (fc F; M*) = * [ d 3 k = 3 ( - (10) 

V ' ^ (27r) 3 7 (k 2 + M* 2 f/ 2 \M* E* F j V ; 

We remark that A(kF, M*) is certainly very small at low densities, and actually it can be still neglected 
up to a baryon density p B ~ 3po (see ref. ]19[). 
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Then in the density range of interest here we can use, at the leading order, a much simpler form of the 
symmetry energy, with transparent 5-meson effects: 



E 



sym 



/ \ 1 fop 1 

{PB) = 6Et + 2 



r i i Ml 

Jp JS \ g* 



Pb 



(11) 



We see that, when the <5 is included, the observed 04 value actually assigns the combination [f p — fs{-^;) 2 ] 
of the {p,5) coupling constants. If fg ^ we have to increase the p-coupling (see Fig.l of ref. p9[). In our 
calculations we use the value 04 = 32MeV. 

In Table I the Set I corresponds to f$ = 0. In the Set II fg is chosen as 2.0/to 2 . Although this value is 
relatively well justified, fl9|| , we stress that aim of this work is just to show the main qualitative new dynamical 
effects of the <5-meson coupling. In order to have the same 04 we must increase the p-coupling constant of a 
factor three, up to f p = 2.8/m 2 . Now the symmetry energy at saturation density is actually built from the 
balance of scalar (attractive) and vector (repulsive) contributions, with the scalar channel becoming weaker 
with increasing baryon density |l9[ . This is indeed the isovector counterpart of the saturation mechanism 
occurring in the isoscalar channel for the symmetric nuclear matter. From such a scheme we get a further 
strong fundamental support for the introduction of the ^-coupling in the symmetry energy evaluation. 




12 3 

Pb/Po 

FIG. 1. Total (kinetic+potential) symmetry energy as a function of the baryon density. 
(NLH + p). Dotted: Hartree (NLH + p + S). Solid: Hartree-Fock (NLHF). 
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FIG. 2. Baryon density variation of the isovector effective coupling when the Fock terms are included. 



Details of the calculation can be found in Refs. [|25| , |l9|] , here we show only the result corresponding to the 
parametrizations (Sets I, II and NLHF) that we will use to investigate the dynamical response. In Fig.|l| 
we report the symmetry energy as a function of the baryon density in the Hartree case (p and p + S) and with 
the Fock terms [(p + 5) case] . When the S "channel" is included the behaviour is stiffer for the relativistic 
mechanism discussed before, see also ref. [[l9| . 

When the Fock terms are evaluated the new "effective" couplings Eqs.(^) naturally acquire a density 
dependence. This is shown in Fig.H for the isovector terms. The decrease of the "effective" p coupling at 
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high density accounts for the slight softening of the symmetry energy seen in Fig]l] (from the dotted to the 
solid line). 



III. LINEAR RESPONSE EQUATIONS 



In this section we study collective oscillations that propagate in cold nuclear matter due to the mean field 
dynamics. In some sense we follow a relativistic extension of the method introduced by Landau to study 
liquid- 3 He 30 32 and recently applied to investigate stable and unstable modes in nuclear matter |20| , |33| , p4j . 
The starting point is the kinetic transport equation Eq.(^). We look for solutions corresponding to small 
oscillations of F(x,p) around the equilibrium value. Therefore we put 



F(x,p)=H(p)+G(x,p) 



(12) 



where H(p) is the Wigner function at equilibrium (see Appendix) and G(x,p) represent its fluctuations. 

For the equilibrium state the Wigner function contains only the isoscalar term and the third component 
of the isovector term. We limit ourselves to study excited states in these channels. Therefore we consider 
isovector density fluctuations with mj- = only, i.e. we do not study processes where a neutron converts into 
a proton or viceversa. In a linear response scheme oscillations in the forementioned channels are decoupled 
from the remaining ones (mj = ±1). 

In the linear approximation, i.e. neglecting terms of second order in G(x,p), the equations for the Wigner 
functions become 



for protons, and 



-d^G {1) {x,p) + (n M - f P b^G {1) (x,p) - M*G {1) (x,p) 
(l- l -A)(f(x)+Mx))H {1) ( P ) , 



-^7 M G(2) (*, p) + (n„ + f p b^G {2) (x, p) - M* G m (x, p) 



(13) 



{l-~A)(f{x)-h{x))H {2) {p) , 



(14) 



for neutrons, where = M — f a ps — fs PS3 and M| = M — f a ps + fs PS3- The quantities jF(x) and T 3 (x) 
are the isoscalar and the isovector components of the self-consistent field: 



T{x) = -8f a G(x) + 8f^G^(x) - 8^ Ps G(x) - 8^ Ps G»(x) 

ops O]^ 

-8-^ Ps G 3 (x) - 8^ Ps G»{x) + 8^ 1 ^G(x) , 
OPS3 oj 3pi dps 



(15) 



A(x) = Sf 5 G 3 (x) + 8f p >y M G$(x) 



PS3 G(x)+8^ 1 ^G(x) 
ops ops 



(16) 



The Hartree approximation is recovered by vanishing all the derivaties inside the quantities jF(x) and 
^(x), except §£|, since still f a = $(p s )/ps ■ 

In order to obtain the equations for the collective oscillations we multiply Eqs. ( |l3| ) and ([lj]) by 7a. 
After performing the traces, and we equate to zero both the real and imaginary parts of the result p0| , pl| . 
Furthermore, by Fourier transforming and integrating over four-momentum, we get the set of equations for 
the scalar and vector fluctuation of each species ( i = 1, 2 for proton, neutron respectively): 



E 

3=1 



S t , 3 + [^~4C^(k))D? 1+ 4C^\k)B 



v 

m ij 



G U) (k) 



4CW(fe)l>TJ + (M_4C7W(fc)) B 



s 



(17) 
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with 



E 

3=1 



(18) 



Aw = /<r + + ^PS ~ M)^PS - (-iy^ P3 



dp 



dp 3 



dpt 



Dl = f u + (-l) i+j f P + ^Ps ~ (-1)#P3 



dps 



dps 



BY...= 



dp. 



-In 



{-i) ~ — Jan 
dps 



B s =2jk 0q -(-iY-iZ 
,. djv n: " 



,df± 



The explicit expressions ol the coefficients (fc), Cjj*' (fc), and C^l(k), together with some of their properties 
can be found in the Appedix. 

We notice, however, that the formalism developed in Hartree-Fock approximation allows to achieve the 
set of equations valid also for an approach to QHD based on a dependence on the scalar density of all 
couplings [p5[ . In particular in our case the coupling functions to the scalar-isoscalar channel depend on all 
the isoscalar and isovector densities and currents (for details, see (§4|), then it is even more general. 

The set of equations developed in the Hartree-Fock approximation recover the ones correponding to the 
usual Hartree approximation (RMF). As already mentioned, it is easly obtained by considering the coupling 
fi to each channel equal to the coupling constant of the corresponding meson. The result is an appreciable 
plainer structure of the equations due to the constant value of all couplings, except f a |36|. The Eqs. ( |l7| ) 
and (ffa) can be reduced to : 



E 



J i,3 



fg-4<7«(fc)K 



(19) 



E 

j'=i 



with Df- - 

6 J 



— + (- 



PS* A 

M* M 



M* 



D 



G»{k)-Uc^\k)D! 3 ] G U) (k) =0 



(20) 



v % = f u + {-iy+>f p . 

The physics effects appear more trasparent and we will follow the Hartree scheme in the next sections 
keeping well in mind that the Fock contributions can be easily included. We expect to have some extra 
contributions in the various interaction channels without qualitative modifications of the physical response. 

The normal collective modes are plane waves, characterized by the wave vector (fc^ = (fc°, 0, 0, | k |)), they 
are determined by solving the set of homogeneous Eqs.( p^ , p0| ). The solutions correspond only to longitudinal 
waves and do not depend on k° and |k| separetly, but only on the ratio 

k° 

Vs = w\- 

The sound velocities are given by values of v s for which the relevant determinant of the set (Eqs.(|l7|, |l8|)) 
vanishes, i.e. the dispersion relations. In correspondence the neutron/proton structure of the eingenvectors 
(normal modes) can be derived. It should be remarked that in asymmetric nuclear matter isoscalar and 
isovector components are mixed in the normal modes. Here this can be argued by the fact that in each of 
the Eqs. ([jl], |2(]) both proton/neutron densities and currents are appearing. 

However we remind that one can still identify isovector-like excitations as the modes where neutrons and 
protons move out of phase, while isoscalar-like modes are characterized by neutrons and protons moving in 
phase (200. 
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IV. THE ROLE OF SCALAR/VECTOR FIELDS IN THE DYNAMICAL RESPONSE 

Before showing numerical results for the dynamical response of asymmetric nuclear matter, in various 
baryon density regions and using the different effective interactions, we would like to analyse in detail the 
structure of the rclativistic linear response theory in order to clearly pin down the role of each meson coupling. 



Isovector Response 

One may expect that once 04 is fixed, the velocity of sound is also fixed |l6[ ] . On the other hand our results 
clearly indicate a different dynamical response with or without the (5-meson channel, for interactions which 
give exactly the same 04 parameter, see Figs.(^|, |]) in the following. In order to get a clear understanding of 
this effect we will consider the case of symmetric nuclear matter in the Hartree scheme, where the dispersion 
relations are assuming a transparent analytical form. 

For symmetric NM the densities, the effective masses and the coefficients C^(k), C^\k) and C^(k) are 
equal for protons and neutrons. Now it is also possible to decouple the collective modes into pure isoscalar 
and isovector oscillations, ^0|. After a straightforward rearragement of Eqs. (|f9|), (p0|), we have for the 
isovector modes: 



5p 3 - 
8ps3 + 



P^_ 8 Cbo(fc) (1 _ 2) 
M* M* V sJ 



Ps 



f P Sp 3 -Sf 5 C (k)Sps3 = 



M* 



8 C(k) f s 5 PS3 + 8f p Co(fc) (1 - v 2 s ) Sp 3 = 



(21) 



We stress that the structure is the same for the isoscalar excitations, of course one has to change the isovector 
fluctuations with the isoscalar ones (Sps, Sps), and the coupling constants of isovector mesons with those 
of the isoscalar mesons . 

Note that in this case to find the zero-sound velocity one has to evaluate determinant of a 2 x 2 matrix 
(and not a 4 x 4), hence the condition for having a solution can be written as 



1 + Np 



/ P (l - vl) - fs^r (l ~ fsA(k F , M*) 



E* 2 



fp 



PS 

M* 



<p(s) = 



(22) 



here Np 



2K F E* 



is the density of states at the Fermi surface and s = v s /vp. To get Eq.(|22|) we have 



used the expression for C(k), Co(fc) and Coo(k) in terms of the Lindhard function ip(s) (see Appendix). The 
quantity A(Kp, M*) is the same integral discussed in Eqs. (|g| JTo| , |TT| ) . 
At this point we can make the following approximation 



Ef 



to evaluate the expression inside the square brackets. Looking at Fig.s (§,§,0) this is a good approximation 
within a 3%. The Eq.(E3) asssumes a quite clear form: 



6E* F 
kp 



sym 



2 Ef 



M* 
771*2 P s 



PB 



if(s) = . 



(23) 



where the potential part of the symmetry energy explicitly appears in the dispersion relations, but joined to 
an important correction term which shows a different f P ,fs structure with respect to that of E^"^, Eqs.(|^, 
|TT| ) . We can easily have interactions with the same a^ value at normal density but with very different isovector 
response. E.g. when we include the S channel we know that we have to increase the f p coupling in order 
to have the same 04, see the discussion of the Eqs.(^, [ll]), but now the "restoring force" (coefficient of the 
Lindhard function in the Eq. (p3|) ) will be strongly reduced. 

Equation (p3|) suggests to define an effective symmetry energy like 



E* = E pot 

sym sym 



2 Ef 



1 - fs Ps 



M* 
Ef 



PB- 



(24) 



which acts as a restoring force for the isovector mode. We can see that once the symmetry energy is fixed its 
effect on the dynamical response depend on the strength of each isovector field. In particular we can easily 
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verify that the fs factor inside bracket in Eq.([24j) is a second order correction and the "leading contribution" 
to the reduction AE* = — E* ym is essentially given by the coupling of the p-meson field. On the 

other hand we know j p5| that once the symmetry energy at saturation density 04 is fixed, the change of 
f p only due to the strenght of fg. We have seen from Table I that /„ can go from 1 fm 2 , if we switch off 
the 5-channel, to f p = 2.8 fm 2 . In terms of the effective symmetry energy this means (if we consider the 
dynamical response at po), AE* ym ~ 4 MeV if fs ~ 2.0 fm 2 . This "softening" of the restoring force easly 
accounts for the decrease of the sound velocity [v s /vp n seen in Fig.^|) for symmetric nuclear matter, a = 0, 
when we pass from NLH — p and NLH — (p + 5). 



Isoscalar Response 

As already remarked, we like to note that for symmetric NM there is a tight analogy between the isoscalar 
and the isovector response in the RMF approach. In the isoscalar degree of freedom the compressibility will 
play the same role of the symmetry energy in the dispersion relation equations. Also in this case we will 
have an important correction term coming from the interplay of the scalar and vector channel. 
(22]) now becomes |36|: 



The Eq. 



l + iV P 
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that can be reduced to the isoscalar equivalent of the Eq. <Y2 
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(25) 



(26) 



where the K^ M is the potential part of the nuclear matter compressibility that in the Hartree scheme has 
the simple structure [^6|(see also Eq.(16) of ref. 
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(27) 



By means of such an analogy, the previous discussion can be extended to isoscalar oscillations with the role 
of E sym now "played" by the compressibility. In this case however one always takes into account both the 
scalar and vector channel in any RMF models. However the coupling costant f u can assume very different 
values depending on the required value for effective masses Mq. This is easy to understand since in the 
RMF limit the saturation binding energy has the simple form 

E/A(0)=E F + f uPB (0)-M N 

where Mm is the bare nucleon mass. So we see that the same saturation values of pb, E/A when decrease 
Mq we have to increase f w . Just to get an idea, we mention that among the most common used RMF 
parametrizations, f u can go from 3.6 fm 2 of NL2 || to 10.2 fm 2 of NL3 |12"| , just decreasing the effective 
mass at saturation from 0.82A/at to 0.6A/jy. In some sense this is obvious because if two EOS have different 
effective masses even if the compressibility is equal the dynamical behaviour is expected to be different. This 
is a very general feature present also in non-relativistic approaches. 

From studies on monopole resonances in finite nuclei with RMF it seems that a higher value of compress- 
ibility is required respect to non-relativistic calculations. Many authors state that this certainly demands 
for a clarification Jl^,|ljl . Even if the monopole resonance is not directly connected to the isoscalar collective 
mode in nuclear matter, our discussion nicely suggests to look at the interplay between effective mass and 
compressibility. For example we can estimate by means of Eq. ©that we can have a shift between the com- 
pressibility and the "effective compressibility" of the order of ~ 100 M eV among different parametrization 
with the same Kmm- Therefore model with K ~ 300 AMeV can reproduce the same frequencies of other 
models with K - 200 AMeV (and a slightly larger M *). 



Landau Parameters 



We would like to briefly discuss the relativistic equations for collective modes in terms of the Landau 
parameters. Interesting features will appear from the comparison to the non-relativistic analogous case. We 
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will focus first on the isovector response, but, as already clearly shown before, the structure of the results 
will be absolutely similar in the isocalar channel. 

The general non-relativistic expression for the isovector modes can be found in Ref. |3l| : 



1 



pa 

F a + 1 c 2 
*° + 1 + 1/3 Ff S 



tp(s) = (28) 



where F£ is the "isovector" combination of the Landau i*o parameters for neutrons and protons Fq = 
pnn _ p™P f that can be expressed in terms of density variations of the chemical potentials: 

Fr = ^N q - 5 qq , , N q = ^1 , q = n, P (29) 

F® are the equivalent for the momentum dependent part of the mean field. In the relativistic approach, for 
symmetric nuclear matter, we get: 

M* 2 1 
FS = F P -F S £*2 1 + fsA (k F ,M*) 

*? - -^TTW ' (30) 

3 P L F 

where Fi — Nf/z (i = p, 5) with Nf = 2N n ._n. Note that the F" contribution comes only from the vector 
coupling. By using the expression Ef"^ Eq.(||), we can write Eq.(28) in the same form of Eqs.(p2], p3| ). The 



result is a similar expression but with the lack of the term in fs inside the brackets in Eq.(|23|). As said, this 
is not the leading term, however around saturation density it amounts to about a 10% of the total correction. 

Moreover turning to the analogy with isoscalar channel the coupling of the a field is now much larger and 
this purely relativistic contribution could be up to a 20%. We underline this point because generally the 
linear response in RMF is discussed calculating the Landau parameters and then using these estimations 
directly into the non-relativistic expression for collective modes Jl6| , [l8| . 

In conclusion from the analysis in terms of the Landau parameters, we can describe the effect of the 
scalar-vector coupling competition previuosly discussed in the following way. The symmetry energy fixes the 
F§, in fact: 

E sym = -^r(l + F a ) (31) 

but in the dynamical response enters also the F^, linked to the momentum dependence of the mean field, 
mostly given by the vector meson coupling. The results are completly analogous in the isoscalar channel, 
with the compressibility given by 

K NM = ^f(l + F*) (32) 

F 

with "isoscalar" combination F* = F™ + Fq P . The relativistic forms of the isoscalar Landau parameters 
are exactly the same as in Eq.(^o|), just substituting the 5,p coupling constants with the a, to ones [ j36[| . 



V. ISOVECTOR COLLECTIVE MODES IN ASYMMETRIC NUCLEAR MATTER 

In this section we discuss results for the isovector collective oscillations which are driven by the symmetry 
energy terms of the nuclear EOS. The aim is mainly to investigate the effect of the scalar-isovector channel. 
This is normally not included in studying the isovector modes and in general the properties of symmetric 
matter in a relativistic approach, while it should be naturally present on the basis of the analysis shown in 
the previous section (and in ref. |19| for equilibrium properties) . Moreover we stress again that Hartree-Fock 
scheme embodies in any case the presence of a scalar-isovector channel, even without the inclusion of the 
5 -meson field |25[ |. 

We will first show results obtained in the Hartree scheme (NLH) including either both the isovector p 
and S mesons or only the p meson. Even if the Hartree approximation has a simpler structure, it contains 
all the physical effects we want to point out. Finally from the complete Hartree-Fock (NLHF) calculations 
we will confirm the dynamical contribution of the scalar isovector channel. 
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For NLH calculations we use the parametrizations of Table I, Set I (p) and Set II (p + S). In the Hartree- 
Fock case the coupling constant fs is adjusted to the value fs(po) = 2.0/m 2 of the NLHF model Eq.(|^). 
We note that this value is smaller that the prediction of recent Dirac-Brueckner-Hartree-Fock calculations 
p8[ , therefore the effects due to the ^-channel presented in the following could be even underestimated in 
the NLHF case. 



Hartree Results 



Let us start by considering isovector-like excitations. In Fig.ga we show the sound velocities in the Hartree 
approximation, as a function of the asymmetry parameter a for different baryon densities. We actually plot 
the sound velocities in units of the neutron Fermi velocities. This is physically convenient: when the ratio is 
approaching 1.00 we can expect that this "zero" sound will not propagate due to the strong coupling to the 
"chaotic" single particle motions ( "Landau damping" ) . This quantity then will also directly give a measure 
of the "robustness" of the collective mode we are considering. 

Dotted lines refer to calculations including (p + 5) mesons, long-dashed lines correspond to the case with 
only the p meson. Calculations are performed at ps = Po and ps = 2 po- We stress that the results of 
the two calculations differ already at zero asymmetry, a = 0. At normal density (po curves), in spite of the 
fact that the symmetry energy coefficient, 04 = E sym (po), is exactly the same in the two cases, significant 
differences are observed in the response of the system. From Fig.0(a) we can expect a reduction of the 
frequency for the bulk isovector dipole mode in stable nuclei when the scalar isovector channel (6 — like) 
is present. Moreover we note that, in the NLH — p case, the exicitation of isovector modes persists up to 
higher asymmetries at saturation density. 




0.2 0.4 

a = (N-Z)/A 

FIG. 3. Isovector-like modes: (a) Ratio of zero sound velocities to the neutron Fermi velocity Vfu as a function of 
the asymmetry parameter a for two values of baryon density. Long dashed line:NLH — p; Dotted line:NLH — (5 + p). 
(b) Corresponding ratios of proton and neutron amplitudes. All lines are labelled with the baryon density, 
po = 0.16/m~ 3 . The full circles in panel (b) represent the trivial behaviour of — (p P /pn) vs. a. 



These are non-trivial features, related to the different way scalar and vector fields are entering in the 
dynamical response of the nuclear system. Such behaviours are therefore present in both collective responses, 
isoscalar and isovector. We have devoted the whole previous section (Sec. IV) to a complete discussion of 
this effect. 

Differences are observed even at ps = 2po, where however also the symmetry energy is different. A larger 
E sym is obtained in the case including the 8 meson (see Fig.|l|) and this leads to a compensentation of the 
effect observed at normal nuclear density. In particular, at higher asymmetries a the collective exicitation 
becomes more robust for NLH — (p + 5). Differences are observed also in the " chemical" structure of the 
mode, represented by the ratio 6p p /Sp n , plotted in Fig.||(b). The ratio of the out of phase n,p oscillations 
is not following the ratio of the n,p densities for a fixed asymmetry, given by the full circles in the figure. 
We systematically see a larger amplitude of the neutron oscillations. The effect is more pronouced when the 
8 (scalar- isovector) channel is present (dotted lines). 
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Hartree-Fock Results 



We have also performed the calculation in the more general case of the Hartree-Fock approximation, 
NLHF, whose formalism has been presented before, Eqs.([l7|, 18). We have fitted the same properties of 



symmetric NM at the saturation density as for Hartree case, NLH. In particular at po the value of the 
isovector coupling is fixed in order get the same symmetry energy (the 04 parameter) of the NLH — (p + 5) 
case. 
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FIG. 4. The same as Fig.3 but for NLHF (solid lines) and NLH - (S + p) (dotted lines). 



In Fig.[| we can see that quite similar results are obtained in Hartree-Fock calculations, with respect to 
the Hartree results including p and 5 mesons, especially at the normal density. This can be understood by 
considering that in Hartree-Fock calculation the effective density dependent couplings associated with the 
isovector channels are tuned in such a way to roughly reproduce, at normal density, the values of the coupling 
constants f p and fs of the Hartree scheme: then not only 0,4 is the same but also its internal structure. Since 
such a tuning can be done only at a given density value, some differences are observed at ps = 2poj due 
to the density dependence of the effective coupling constants of the NLHF scheme. Therefore the greater 
value of sound velocity in the Hartree case can be linked to the greater value of E symi see Fig. Ill 



Disappearance of the Isovector Modes 



For asymmetric matter we have found that, in all the calculation schemes, with increasing baryon density 
the isovector modes disappear: we call such densities p c g oss . E.g. from Figs.|](b),f|(b) we see that the ratio 
Spp/Sp n tends very quickly to zero with increasing baryon density, almost for all asymmetries. Around this 
transition density we expect to have an almost pure neutron wave propagation of the sound. Here we show 
the results of the NLH + p case, see Figs.|| and [| but the effect is clearly present in all the models. 

For symmetric matter we have a real crossing of the two phase velocities, isoscalar and isovector, as shown 
in Fig.g(a). Above p c g oss the isoscalar mode is the most robust. 
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FIG. 5. Sound phase velocities of the propagating collective mode vs. the baryon density (NLH + p case). Crosses: 
isovector-like. Open circles: isoscalar-like. (a): symmetric matter, (b): asymmetric matter, a = 0.1. (c): asymmetric 
matter, a = 0.5. 
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FIG. 6. Ratio of protons and neutron amplitudes in the propagating mode, for different asymmetries, as a function 
of the baryon density around the pg'° as . Crosses: a = 0.1, Fig. 5b. Open circles: a = 0.5, Fig. 5c. 



For asymmetric matter we observe a transition in the structure of the propagating normal mode, from 
isovector-like to isoscalar-like, Fig.H(b,c). Similar effects have been seen in a non-relativistic picture pofl . 

For a given asymmetry a the value of p c g oss is different for the three models considered, as can be argued 
by the behaviour of 8p p /5p n at 2,00 in Fig s -||(b), @(b). E.g. for a = 0.1 NLHF has the lower value 
[p c g oss — 2.4p ), while NLH — p has the higher one (pq oss — 3.0po)- This is again related to the reduction 
of the isovector restoring force when the scalar-isovector channel (5-like) is present, see Sec. IV. 

From Fig.g we see that the proton component of the propagating sound is quite small in a relatively wide 
region around the "transition" baryon density, a feature becoming more relevant with increasing asymmetry, 
see the open circle line. This is quite interesting since it could open the possibility of an experimental 
observation of the neutron wave effect. 



VI. ISOSCALAR COLLECTIVE MODES IN ASYMMETRIC NUCLEAR MATTER 

So far we have focussed our discussion on the isovector-like response of the asymmetric nuclear matter. 
However it is well known that in asymmetric nuclear matter can exist also isoscalar-like modes, see [ p0||3^ ] 
and rcf.s therein. 



Exotic High Baryon density modes 

From the previous analysis we have seen the isoscalar-like excitations to become dominant at high baryon 
density, above the p c £ oss introduced before. 
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Some results are shown in Fig]?]. It should be noticed that the frequency of the isoscalar-like modes is 
essentially related to the compressibility of the system at the considered density. In Fig.^(a) we display 
the sound velocity obtained in Hartree and Hartree-Fock calculations at ps = 3.5 po, as a function of the 
asymmetry a. The differences observed among calculations performed within the Hartree or Hartree-Fock 
scheme are due to a different behaviour of the associated equation of state at high density. 




a = (N-Z)/A 

FIG. 7. The same of Fig. 3 for isoscalar-like modes, at ps = 3.5po- Solid line: NLHF. Long Dashed line: NLH — p. 
Dotted line: NLH ~ (p + 5). The full circles in panel (b) represent the behaviour of p p /p n vs. a. 

At a — the two Hartree models have exactly the same isoscalar mean fields, but for asymmetric nu- 
clear matter the different behaviour of the symmetry energy leads to a different compressibility The case 
NLH — (p + S) which has the stiffer E sym (resulting in a greater incompressibility for a > 0) with respect 
to NLH — (p) shows also a greater increase of v s /vfu with density. Instead, NLHF even if it has the 
same compressibility Knm at saturation density, shows a different v s . This should be due to the density 
dependence of the coupling function arising from exchange terms which leads to different values of Knm out 
of pa (even for a = 0). 

Some differences are observed also in the chemical composition of the mode (Fig.|^(b)). The black spots 
show the behaviour of p v j p n vs. a. Note the pure neutron wave structure of the propagating sound, since 
the oscillations of protons appear strongly damped (Sp p /Sp n <C p p /p n ); unfortunately this is an effect not 
experimentally accessible (at present), see also the discussion at the end of the previous Section. 

Before closing this discussion we have to remark that the isoscalar-like modes at high baryon density are 
vanishing if the nuclear EOS becomes softer. This is indeed the results of two recent models, ref.s p8|p9[ , 
where the nuclear compressibility is decreasing at high baryon density for a reduction of the isoscalar vector 
channel contribution. In fl8|| this is due to self-interacting high order terms for the ui meson, while in ref. 
p9[ to a reduced f u coupling with increasing baryon density. 

Finally we note that all causality violation problems (superluminal sound velocities) observed in the non 
relativistic results at high baryon density, see Jl6| and Fig. 3c in ref. |2(J], are completely absent in the 
relativistic approach, see the high density trends in FigM 



Isospin Distillation in Dilute Matter 



We have also investigated the response of the system in the region of spinodal instability associated with 
the liquid-gas phase transition, which occurs at low densities. It is known that in this region an isoscalar 
unstable mode can be found, with immaginary sound velocity, that gives rise to an exponential growth of the 
fluctuations. The latter can represent a dynamical mechanism for the multi-fragmentation process observed 
in heavy-ion collisions. We have found this kind of solution in the present approach. In Fig.^ we show the 
ratio 5pp/Spn as function of the initial asymmetry for such a collective mode. For all the interactions this 
ratio is different from the corresponding p p /p n of the initial asymmetry a. This is exactly the chemical effect 
associated with the new instabilities in dilute asymmetric matter [H,p7|. 
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FIG. 8. Isocalar-like unstable modes at ps = 0.4 po ■ Imaginary sound velocity (a) and ratio of proton and neutron 
amplitudes (b) as a function of the asymmetry a. Solid line: NLHF; Dotted line: NLH — (p + 5); Long Dashed 
line: NLH — p. The full circles in panel (b) represent the behaviour of 5p p /8p„ vs. a. 

In particular it is found that, when isoscalar-like modes become unstable, the ratio Sp p /Sp n becomes 
larger that the ratio p p /p n (at variance with the stable modes at high densities, see Fig]?]). Hence proton 
oscillations are relatively larger than neutron oscillations leading to a more symmetric liquid phase and to a 
more neutron rich gas phase, during the disassembly of the system. This is the so-called isospin distillation 
effect in fragmentation, and signatures of this effect could be searched by looking at the ratio N/Z of 
fragments produced in dissipative heavy ion collisions pO| , ^l| . 

We note here that in dilute asymmetric NM we can distinguish two regions of instability, mechanical 
(cluster formation) and chemical (component separation) . There is however no discontinuity in the structure 
of the unstable modes which are developing. For all realistic effective nuclear interactions (relativistic and 
non) the nature of the unstable normal modes at low densities is always isoscalar-like, i.e. with neutrons 
and protons oscillating in phase, although with a distillation effect discussed before, see ref. |37 for a fully 
detailed study of this important property of asymmetric nuclear matter. 

In Fig.|| we observe that Hartree results (with and the without 8 meson) are very similar and, indeed, at low 
density the symmetry energy behaviour is nearly the same in the two cases. On the other hand, differences 
are observed with respect to the Hartree-Fock case. In fact, the NLHF symmetry energy presents a softer 
behaviour (around ps = 0.4po), that leads to a smaller distillation effect. We remark that the equality 
between NLH with and without 5 is in agreement with the analysis in terms of the generalized Landau 
parameters associated to normal modes developed in Ref. [fiofj . We can conclude that there are essentially 
no effects of the scalar isovector channel on isospin distillation in the spinodal decomposition. 



VII. CONCLUSIONS AND OUTLOOK 



We have developed a linear response theory starting from a relativistic kinetic equations deduced within a 
Quantum-Hadro-Dynamics effective field picture of the hadronic phase of nuclear matter. In the asymmetric 
case we consider as the main dynamical degrees of freedom the nucleon fields coupled to the isoscalar, scalar 
a and u>, and to the isovector, scalar S and vector p, mesons. 

Using the Landau procedure we derive the dispersion relations which give the sound phase velocity and the 
internal structure of the normal collective modes, stable and unstable. We have focussed our attention on 
the effect of the isovector mesons on the collective response of asymmetric (neutron- rich) matter. In order 
to better understand the dynamical role of the different mesons, the results are obtained in the Hartree 
approximation, which has a simpler and more trasparent form. The contribution of Fock terms is also 
discussed. 

We have singled out some qualitative new effects of the (S-meson-like channel on the dynamical response 
of ANM . Essentially, our investigation indicates that even if the symmetry energy is fixed, the dynamical 
response is affected by its internal structure, i.e. the presence or not of an isovector-scalar field. This is 
implemented by the explicit introduction of an effective 5-meson and/or by the Fock term contributions. 
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The interest comes from the fact that both mechanisms are absent in the usual relativistic RPA calculations 
for finite nuclei. 

It is important to stress that the same interplay between scalar (er-meson) and vector (w-meson) contri- 
butions can be seen in the dynamical isoscalar response. In general we clearly show a close analogy in the 
structure of the linear response equations: 

• Same form of the dispersion relations, cfr. Eq.(p2|) and Eq.(|25|). 



Parallel role of E^ m and K^ M in the determination of the restoring force, Eqs.(J23h and (27) 



Parallel structure of the corrections due to the scalar- vector meson competition, Eqs.(23) and ( 



This appears to be a beautiful "mirror" structure of the relativistic approach that strongly support the 
introduction of a 5-meson-like coupling in the isovector channel. We like to remind that the same "mirror" 
structure of the relativistic picture has been recently stressed in ref. for equilibrium properties, saturation 
binding and symmetry energy, the a\ and 04 parameters of the Weiszaecker mass formula. 

The relativistic dispersion relations have been compared with the non-relativistic ones of the Landau 
Fermi Liquid theory expressed in terms of Landau parameters evaluated in a relativistic scheme Jig ] . The 
corrections appear to be not negligible, particularly for the isoscalar response. 

From the numerical results on the collective response of ANM some general features are qualitatively 
present in all the effective interactions in the isovector channel: 

• In asymmetric matter we have a mixing of pure isoscalar and pure isovector oscillations which leads to 
chemical effect on the structure of the propagating collective mode: the ratio of the neutron/proton 
density oscillations Sp n /Sp p is different from the initial p n /p p of the matter at equilibrium. However 
we can still classify the nature of the excited collective motions as isoscalar — like (when neutron and 
protons are oscillating in phase) and isovector — like (ou t of phase). To note that similar effects can 
be obtained also using non-relativistic effective forces |20| , |37|| . 

• For a given asymmetry the isovector-like mode is the most robust at low baryon density, always showing 
a larger neutron component in the oscillations. With increasing baryon density we observe a smooth 
transition, at a p c g oss ~ 2 — 3poi to an isoscalar-like branch, still with a dominant Sp n . In the region of 
the transition we predict a propagation of almost pure neutron waves. For relatively large asymmetries 
(a = fej = 0.5, N = 2>Z) this behaviour is present in a wide interval of densities around p c £ oss . All 
that seems to suggest the possibility of an experimental observation of related effects in intermediate 
energy heavy ion collisions with exotic beams. If the compressibility of nuclear matter is decreasing at 
high baryon density also these exotic isoscalar-like mode will disappear. This could be a nice signature 
of the softening of nuclear EOS at high densities. 

• The isoscalar-like motions become unstable at sub-saturation densities still with a strong chemical 
effect, now in the opposite direction with respect to the one discussed before, present in the stable high 
density modes. Now the unstable oscillation is more proton-rich, eventually leading to the formation 
of more symmetric clusters vs. a very neutron-rich gas phase. This is the neutron distillation effect 
J^ p0]p^j3^7| , [f0| . [4l| | , a new important feature of the liquid-gas phase transition in asymmetric nuclear 
systems. 
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VIII. APPENDIX 

The Wigner matrix H(p) for matter at equilibrium saturated in spin has the following form: 

^From the kinetic equation one obtain the relation for nuclear matter at equilibrium between the scalar and 
vector parts: 



1G 



where the zero component of the vector part is proportional to the Fermi Dirac distribution function: 



1 1 



*(*) 



where i = n,p and E* t = (k 2 + M* 2 ) 1 ' 2 . 

The coefficients C«(fc), C%\k) and C^(fc) introduced in Sec. II are given by the integrals 
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= j*p 



H'^(p)Kp»p; 

p* p l kP 



(Ale) 



where H'A(p) — d^H^(p). The index i specifies the kind of nuclcon: i = 1 for protons and i = 2 for 

neutrons. The frequency fc° includes an imaginary part ie with e positive infinitesimal. 
By using the definitions (Al) it can be easily checked that 

fc*cf(fc) = 0, ^cf>(k) = -k/-f, 



Y J C^(k) = M*C^(k)-\p Si 



(A2) 



In order to be more specific we choose the z axis in the direction of the wave vector k. As a consequence, 
the following coefficients identically vanish: 

Cf (fc), Cf (fc), C«(fc), C«(fc), and C«(fc) 

for I ^ m (I and m are space indices ). In addition, for symmetry reasons, 

cffCO = C22 w • 

The integrals in Eqs. (Al) can be evaluated analytically. They give 



cW=+JiV^( Si ) , 



C11 CO = 4 A»« - g ~j^~P Bi + 8EpS Sl ~ ^ ^ 



(A3a) 

(A3b) 
(A3c) 
(A3d) 



where vp t is the Fermi velocity, Sj = k°/(vp i |k|), iVj are the density of states at Fermi surface and 



tp(si) = 1 - -^-Zn 



7TS;0(1 - Si) 



is the Lindhard function. The remaining coefficients C^(fc), C^(fc) and C^(fc) can be evaluated by means 
of the relations (A2). 
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